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The magnetic phase diagram of PrOs4Sbi2 has been investigated by specific heat measurements 
between 8 and 32 T. A new Schottky anomaly due to excitations between two lowest crystalline- 
electric-field (CEF) singlets, has been found for both H \\ (100) and H | (110) above the field where 
the field-induced ordered phase (F10P) is suppressed. The constructed H — T phase diagram shows 
weak magnetic anisotropy and implies a crossing of the two CEF levels at about 8 - 9 T for both field 
directions. These results provide an unambiguous evidence for the Fi singlet being the CEF ground 
state and suggest the level crossing (involving lowest CEF levels) as the driving mechanism of FIOP. 



PACS numbers: 71.27.+a, 74.70.Tx, 75.40.Cx 

A recent discovery of superconductivity 1] in 
PrOs4St>i2 is an important milestone in the field of 
strongly correlated electrons for a number reasons. 
Firstly, a large discontinuity in the specific heat at 
the superconducting transition temperature T c un- 
ambiguously establishes that Pr-based materials can 
support heavy fermion (hf) states. Secondly, proposed 
crystalline electric field (CEF) schemes suggest non- 
magnetic ground states and exclude a conventional 
Kondo effect, believed to be the source of hf behavior 
in Ce- and U-based metals. The only other serious 
microscopic mechanism considered as the origin of 
the hf behavior is a controversial quadrupolar Kondo 
effect. |2j Finally, two superconducting transitions 
observed[3j in this compound imply an anisotropic order 
parameter. Therefore, PrOs4Sbi2 is a candidate for a 
novel mechanism of superconducting pairing that is due 
to neither phonons nor magnetic interactions. 

Correct interpretations of these new exciting phenom- 
ena hinge crucially on the proper accounting for the CEF 
configuration of Pr. The consensus is that the non- 
magnetic (non-Kramers) doublet is the ground state, 
and T$ is the first excited triplet. However, available 
experimental data do not rule out the possibility of the 
Ti singlet as the ground state. Among the strongest ar- 
guments for both configurations were those provided by 
specific heat in zero and small magnetic fields. In partic- 
ular, the zero field Schottky-like anomaly at 3.1 K can be 
related to the T$ - r$ model, assuming these two levels 
are split by 7.5 K, or Ti - r$ model with the splitting of 
8.4 K. Inelastic neutron scattering measurements |4| yield 
the lowest-energy peak at 8 K. Entropy changes in small 
fields and as a function of temperature were also used to 
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argue for either of the models Q, H H ||. The difficulty 
with interpreting these low temperature, low field results 
is related to strong hybridization of 4f and conduction 
electrons, as evident from the large electronic specific 
heat coefficient and the size of the discontinuity in the 
specific heat (C) at T c . Therefore, we have performed 
specific heat measurements in high magnetic fields, in 
which a coupling between conduction and f electrons is 
expected to be suppressed, to get insight into the CEF 
configuration of Pr in PrOs4Sbi2. 

The crystal used in our investigation was from the 
same batch of high quality samples for which a host of 
other measurements were reported. 0,0] The specific heat 
study to 14 T was performed in a superconducting mag- 
net. The measurements above 14 T were carried out in 
a resistive Bitter magnet. The calorimetry employed the 
relaxation method, in which a small sample of about 1.5 
mg was used to optimize measuring conditions. Figs. 1 
and 2 show specific heat in fields ranging from 10 to 32 T. 
The field was applied along the nominal (100) direction. 
Our estimate of the probable error in the sample align- 
ment is about 5°. All the data presented are after sub- 
tracting the phonon background approximated by a De- 
bye temperature (8_d) of 165 K (after Vollmer et al.) 3]). 
This value of 0£> is quite controversial. In fact, Bauer 
et al.0, Aoki et al. and Maple et al.0 suggested 
Od to be 304, 320, and 259 K, respectively. Changing 
Od to 320 K does not alter the main conclusions of our 
study, but has some quantitative impact on the results, 
as discussed further. A relatively high low-temperature 
limit was chosen to avoid complications associated with 
a nuclear contribution of Pr. This nuclear contribution is 
strongly enhanced by coupling with orbital moments of 
f-electrons.0 The specific heat, at temperatures where 
the nuclear degrees of freedom dominate, is difficult to 
measure by a conventional relaxation method because of 
an additional time scale (s) entering the experiment, the 
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FIG. 1: Specific heat of PrOs 4 Sbi 2 in 10 and 12 T in the vicin- 
ity of FIOP transition. The magnetic field was applied along 
the (100) direction. Phonon contribution was subtracted as- 
suming Od =165 K. 



nuclear spin-lattice relaxation time T\ . Q Strongly non- 
exponential temperature decays observed at the lowest 
temperatures (e.g., below 0.5 K in the field of 10 T and 
below 1.5 K in the field of 32 T) indicate the importance 
of nuclear degrees of freedom and cannot be analyzed us- 
ing the so called r 2 correction. Therefore, these lowest 
temperature points carry large uncertainty and are not 
analyzed in detail. 

Figure 1 shows the low temperature specific heat for 
10 and 12 T. Pronounced peaks in the vicinity of 1 K cor- 
respond to a transition to a field induced ordered phase 
(FIOP), which was first observed by Aoki et al.Q for 
fields larger than 4 T also applied along the (100) direc- 
tion. The highest field used in this latter study was 8 T. 
This long range order was confirmed by other investiga- 
tions, including the specific heat of Vollmer et al.|3j] and 
magnetization study of Tayama et al. 9] It has been pro- 
posed that the order parameter corresponding to FIOP is 
antiferro-quadrupolar type. This type of order is consis- 
tent with large anomalies in the specific heat and a very 
small value of the ordered (antiferromagnetic) moment 
(about 0.025 fi B at 0.25 K) 10]. Combining our results 
with those of Aoki et al. Q and Vollmer et al. shows 
that the transition temperature (T n ) identified as a posi- 
tion of the peak in the specific heat (C) reaches a max- 
imum value somewhere near 9 T. Similarly, C(T n ) has 
a non-monotonic field dependence, attaining the highest 
value between 8 and 10 T. Notice a large reduction of 
C{T n ) and small but clearly discernible decrease of T n 
between 10 and 12 T. It is interesting to note a relatively 
large discrepancy for C at T n between all three reports for 
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FIG. 2: Specific heat of PrOs4Sbi2 in magnetic fiefds 13, 13.5, 
20, and 32 T applied parallel to (100), after phonon subtrac- 
tion. Notice the appearance of a shoulder at about 1.2-1.3 K 
at 13 T and the disappearance of the FIOP transition at 13.5 
T. The low temperature tails are due to nuclear contribution 
of Pr. 



overlapping fields. This discrepancy is probably due to 
the anisotropic response of the specific heat in magnetic 
fields near T n , discussed below, and some misalignment 
of the sample with respect to the (100) direction. An 
increase of the field from 12 to 13 T (Fig. 2) results in 
a small reduction of T n and large suppression of C{T n ). 
Moreover, a shoulder appears on the high temperature 
side of the FIOP anomaly. The specific heat value at 
this shoulder is about 3400 mJ/K mol. This shoulder 
evolves into a broad maximum for a slightly higher field 
of 13.5 T. Again, the value of C at the maximum is be- 
tween 3400 and 3500 mJ/K mol. In this field we no longer 
detect a signature corresponding to FIOP. Thus, our re- 
sults strongly imply a disappearance of FIOP before T n 
reaches 0. A similarly broad peak exists at all fields stud- 
ied up to at least 32 T. The temperature of this broad 
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FIG. 3: Magnetic field phase diagram of PrOs4Sbi2 for the 
field parallel to (100). Filled squares represent the FIOP tran- 
sition. Open squares correspond to the Schottky anomaly. 
The inset shows a model calculation of the H versus T phase 
diagram assuming the singlet as the ground state. The solid 
fine represents the FIOP boundary; the dashed line corre- 
sponds to a maximum in C (T m = 0.417(5; 8 is the energy 
separation between the lowest CEF levels.) The model Hamil- 
tonian and CEF parameters are identical to those of ref. 10. 



maximum increases with the field strength (Fig. 2). 

The magnitude of the maximum, in fields of 13 T and 
larger, is essentially field independent and ranges from 
3300 to 3500 mJ/K mol. These values are within about 
10% of the maximum value for a Schottky anomaly of 
a two level system with identical degeneracies, ^lj The 
uncertainty of our specific heat measurements in these 
fields (and at temperatures where nuclear contribution is 
small) is about 10%. In addition, there is the aforemen- 
tioned uncertainty associated with a subtraction of the 
phonon term. Increasing the Debye temperature from 
165 K, used in our subtraction, to the other extremal 
value proposed, 320 K, raises the estimate of the elec- 
tronic part of C by about 290 mJ/K mol at 3.5 K. Thus, 
the extracted values at the maximum are well within 
the realistic error bar of the theoretical 3650 mJ /K mol. 
Since the 32 T field is large enough to significantly split 
any degenerate levels, the observed Schottky anomaly 
has to represent excitations between two singlets. The 
temperature of the maximum, T m , is related to the en- 
ergy separation of the two levels 6, T m — 0.4175. An 
extrapolation of T m to T = (Fig. 3) determines the 
field at which the two levels cross. Our best estimate of 
this crossing field, using the lowest fields, is between 8 
and 9 T. 



These results can be used to infer new information 
pertaining to a plausible crystal field configuration of 
Pr. Pr can be modeled by the following single-site 
Hamiltonian[To| . 

H = H C kf - SjMbJH - J(J')3 - Qi{Oi)Oi, (1) 

i 

where Hqef, J and Oj represent the CEF Hamiltonian 
for the cubic Th symmetry, the total angular momentum 
and the i-th quadrupole moment (see Ref. 01 f° r detailed 
definition). Note that the conventionally used Wcef0 
with the Oh symmetry |l2j| is only approximate because 
it ignores the the tetrahedral symmetry of Sb ions.|l3|. 
This latter CEF Hamiltonian describes the Zeeman ef- 
fect roughly but cannot account for the neutron diffrac- 
tion results in the FIOP region dOj. In particular, the 
lower Th symmetry is needed to explain the appearance 
of a small antiferromagnetic moment along the (010) di- 
rection, and not along (100), when the field is applied 
parallel to the (001) direction. (Directions (100) and 
(010) are not equivalent). Using the CEF parameters 
proposed by Kohgip^f for the Ti - T 5 CEF configura- 
tion, T m (with Qi = 0) and the O vz -tyve quadrupolar 
ordering temperature T n are calculated [lj]- Note that 
this model consistently reproduces the inelastic neutron 
datay| , characteristic maxima appearing in the magnetic 
susceptibility and the zero-field specific heat ply, 0- 
The results of our calculations are presented in the inset 
to Fig. 3. 

The measured phase diagram (Fig. 3) and the the- 
oretical one expected for the T\ - model (inset to 
Fig. 3) are strikingly similar. In particular, in both 
diagrams, the crossing field is very close to the one at 
which the transition temperature of FIOP becomes max- 
imum. Some overestimation of T n can be attributed to 
the mean-field nature of the calculations. On the other 
hand, there is less defined correlation between these two 
fields in the alternative - model. Instead, it has 
been argued for this latter model that FIOP is related to 
the crossing of the upper T% and lower T$ levels. 3] Thus, 
the observed correlation between the two characteristic 
fields constitutes our first argument for the T± singlet 
being the lowest CEF level. 

More conclusive arguments regarding the CEF config- 
uration can be obtained from the study of the anisotropy 
of the Zeeman effect. The Zeeman effect is expected to 
be strongly anisotropic for the T3 doublet and almost 
isotropic for the Ti singlet, for the range of CEF pa- 
rameters that are consistent with the susceptibility 0. In 
particular, there is no crossing involving the lowest CEF 
levels when the field is applied along the (110) direction 
in the T3 - r§ model. We emphasize that this predic- 
tion is quite general; it does not depend on whether one 
uses the the CEF Hamiltonians and parametrs proposed 
by Bauer et al.p| or Kohgi et al.[ig. Fi gure 4 shows 
the measured high field phase diagram for this field di- 
rection. For this direction we observe a decrease of T n 
values with respect to the (100) direction for the cor- 
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FIG. 4: Magnetic field phase diagram of PrOs4Sbi2 for H 
parallel to (110). For the definition of symbols see Fig. 3. 



responding fields, consistent with the previous magneti- 
zation measurements 9J. On the other hand, within the 
uncertainty of the measurement, there is no change in 
the position of the Schottky anomaly at 13 and 14 T, as 
expected for the Ti CEF ground state and inconsistent 
with the T3 scenario. Moreover, for the (110) orientation 
we clearly observe the Schottky anomaly already at 12 
T. This lower field limit for the Schottky maximum is 
probably due to competition between the two types of 
anomalies, discussed below, and lower values of T n for 
the (110) direction. A straight line fit for the three T m 
points results in the crossing field value of 8.5 ± 0.5 T. 
This value agrees with our estimate for the (100) direc- 
tion. 

Existence of the crossing field for the (110) direction 
provides an unambiguous evidence for the Ti - model, 
in fact the strongest evidence to date. We stress that a 
small misalignment of the sample with respect to the field 
in either of the measurements cannot explain essentially 
identical crossing fields for both directions. In fact, the 
measured difference in T n values for (100) and (110) di- 
rections provides an additional check of the alignment. 



Similarly to the (100) direction, there seems to be a close 
correlation between the crossing field and the field corre- 
sponding to T n maximum. 

Figures 3 and 4 imply a strong competition between 
the field-induced order and Schottky peak. The FIOP 
transition in the specific heat abruptly disappears before 
T n reaches zero. Precise magnetization^ and de Haas 
van Alphen measurements^, on the other hand, were 
able to map T n as a function of the magnetic field all 
the way to T n w 0. This apparent contradiction can 
be explained by a very small entropy available for the 
FIOP transition above 13 and 12 T for fields parallel 
to the (100) and (110) directions, respectively. Specific 
heat, being a bulk measurement, can be less sensitive 
than magnetization techniques in this situation. A strong 
competition is to be expected in the Ti - r§ scenario. The 
ground state pseudo-doublet formed at the level crossing 
carries both magnetic and quadrupolar moments. Since 
a quadrupolar moment operator does not commute with 
a dipolar one, the quadrupolar interactions leading to 
FIOP compete with the magnetic Zeeman effect. 

Summarizing, wc find an evidence for the level cross- 
ing occurring around 8-9 T. This field value is very close 
to the field at which both T n and C(T n ) become max- 
imum, suggesting that the crossing of the lowest CEF 
levels is the driving mechanism of FIOP. In fields above 
10 T there is a strong competition between FIOP and the 
Schottky peak. Our high field specific heat results pro- 
vide compelling arguments for the Ti singlet being the 
CEF ground state of Pr in this system. This results call 
for a new model of electronic mass enhancement, instead 
of the ordinary quadrupolar Kondo effect, which has been 
discussed mainly in the context of the T3 non-Kramers 
doublet ground state. 
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